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Abstract 

I will review recent advances in the field of blazars, highlighting the contribution of Swift. Together with other operating satellites 
(most notably Fermi, but also AGILE, WISE, Planck) and ground based facilities such as Cherenkov telescopes. Swift was (and 
is) crucial for improving our understanding of blazars. The main advances in the blazar field made possible by Swift includes the 
opening of the time domain investigation, since there are several sources with hundreds of simultaneous optical, UV and X-ray 
data taken at different times; the possibility to measure the black hole mass in very powerful blazars, that show clear signs of 
accretion disk emission; the possibility to classify blazar candidates, through X-ray observations; the finding of the most powerful 
and distant blazars, emitting strongly in the hard X-ray band accessible to SwiftfBAT. All these improvements had and have a great 
impact on our understanding on how relativistic jets are formed and emit, on their power, and on how the heavy black holes in 
these systems first formed and grew. 

Keywords: Active Galactic Nuclei; Black hole physics; BL Lac objects; Radiation processes: non-thermal; y-rays; Accretion 
disks. 


1. Introduction 

After the launch of the Compton Gamma Ray Obser¬ 
vatory, CGRO, we discovered that blazars (i.e. quasars 
with jets pointing nearly at us) were the most impor¬ 
tant class of persistent y-ray emitters. This was some¬ 
what unexpected, despite the fact that the previous y- 
ray satellite COSB already detected 3C 273 as a y-ray 
source [!95|, It was not completely clear that the 
production of y-rays was associated to relativistic jets, 
even if all the necessary ingredients were known since 
the early seventies: superluminal motion and the pres¬ 
ence of relativistic electrons in the source, producing 
synchrotron and self-Compton radiation (the external 
Compton idea was yet to come). 

The discovery of the the strong y-ray emission of 3C 
279 in 1991 by CGRO was soon followed by the reali¬ 
sation that blazars are y-ray emitters as a class, and this 
triggered a frantic phase of theoretical developments 
[hi 0, H, 0,0- At the same time, the (few) 


multi-wavelengths campaigns showed coordinated 
variability of the flux at different frequencies, and this 
made the jet paradigm to shift from a multi-zone jet, 
producing the highest frequencies (y-rays) in the inner¬ 
most regions and the IR-optical further out [63], [59], 
I3~3h, to the much simpler “one-zone” jet in which most 
of the emission was produced in a single region, i.e. 
the same electron population producing the synchrotron 
was also responsible for the high energy flux (but not 
the radio, due to the synchrotron self-absorption). This 
required a strong effort, especially on the observational 
side, because it was not easy to organise multiwave¬ 
length campaigns joining space and ground observato¬ 
ries. 

When Swift was launched what was a dream became 
routine: simultaneous optical, UV, and X-rays observa¬ 
tions became easily accessible and flexible planning al¬ 
lowed to use Target of Opportunity observations to fol¬ 
low extraordinary events. Then, when Fermi joined in, 
we could have a really complete view of the behaviour 
of blazars, and not only of the 3 or 4 brightest ones, but 
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of hundreds of them. 

The Swift and Fermi satellite, together with ground 
based facility like the Cherenkov telescopes, made a 
quantum jump in our knowledge if the physics of 
blazars, and led the way to use them not only to un¬ 
derstand the high energy physical processes that char¬ 
acterise their emission, but to use blazars as a probe of 
the far Universe. 

What follows is a partial view of the recent advances 
in blazar science allowed by Swift. 


2. Multi-wavelength campaigns 


Both planned observations together with other in¬ 
struments and target of opportunity (ToO) observations 
(performed after even a very short notice) have secured 
the optical-UV and X-ray observations of hundreds of 
blazars. Fig. |T|shows the observed spectral energy dis¬ 
tribution (SED) of blazars together with the observing 
band of Swift and Fermi/ LAT. Swift/ UVOT and XRT 
cover the peak of the synchrotron emission of low power 
line-less BL Lacs, the best candidate to be TeV emit¬ 
ters, while Swift/BAT can be more effective in observing 
and even discover powerful flat spectrum radio quasars 
(FSRQ) with broad emission lines at high redshift. For 
this class of objects, Swift/UVOT can observe the ther¬ 
mal emission produced by their accretion disk. 

This thermal component may be elusive for inter¬ 
mediate redshift and intermediate power FSRQs, be¬ 
cause in these objects the beamed non-thermal flux 
can hide the thermal continuum. One example 3C 
454.3 1 107 1. [14], j75], whose thermal emission was 
revealed through the optical and UV monitoring involv¬ 
ing Swift/WOT |7(|. Another example is B3 1633+382 
(z = 1.814), discovered as a y-ray source by CGRO and 
well monitored by Swift, Fermi and AGILE [[77 ]. As Fig. 
[2] shows. Swift was instrumental to reveal the contribu¬ 
tion of the accretion disk. 

Fig. [3] shows data from Swift and NuSTAR of the 
high-redshift blazars PKS 2149-306 (z = 2.345) 
Together, Swift/XKV and NuSTAR cover the 0.3-70 keV 
band. The bottompanel of Fig. QJshows the X-ray SED, 
with the two observations of Swift+NuSTAR together 
with other archival observations. It can be seen that 
Swift is crucial to describe the behaviour of the X-ray 
spectrum, that does not change at low energies, while it 
becomes harder when brighter above ~ 4 keV (13 keV 
rest frame). The addition of the FermifLAT makes clear 
that the hard X-ray behaviour corresponds to a shift in 
the high energy peak frequency, that becomes smaller in 
the (slightly) lower state. In this case we have a behav¬ 
ior opposite to the blazar sequence (see Fig. []}. 



Figure 1. The “blazar sequence”: blazars have SEDs that changes 
according to the bolometric observed jet luminosity. Low powerful 
lineless BL Lacs are “blue”: their synchrotron and Compton hump 
peaks at high frequencies, and the corresponding luminosities are 
about equal. Powerful flat spectrum radio quasars (with broad emis¬ 
sion lines) are redder, and the Compton hump dominates. This has 
been explained as due to radiative cooling: electrons in more powerful 
sources suffer more severe losses, and this limits their typical energies 
to values smaller than the one in low powerful BL Lacs, in which the 
cooling is less severe t35ll . In this respect, the presence or absence of 
the broad emission lines can play a crucial role, since they can largely 
enhance the inverse Compton emission and the corresponding radia¬ 
tive cooling. The indicated yellow vertical stripes correspond to the 
observing bands of Swift and Fermi/LACY. Adapted from f2^1 . (2&\ . 
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Figure 2. The Swift! UVOT SED of the blazar B3 1633+382 (alias 
4C38.41), at different epochs. The flux is the sum of the steep tail of 
the synchrotron jet emission and the thermal component produced by 
the accretion disk. From (H. 
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3. Time domain 

The accumulation of data during the life of Swift im¬ 
plies that a source has the chance to be observed several 
times, with all the three Swift instruments. The most fa¬ 
mous blazars (PKS 2155-304, Mkn 421, Mkn 510, 3C 
454.3) have been observed hundreds of times. All data 
are public, and the Swift archive is a resource for the 
years to come, still to be fully exploited. As in many 
other branches of science, the amount of data is becom¬ 
ing too large to be analysed by humans (or, at least, by 
single humans), and as automatic tools have been de¬ 
veloped for Gamma Ray Bursts, there has been the ini¬ 
tiative of the ASI Science Data Center (ASDC) to offer 
tools to build the SED of all sources (not only blazars) 
with the option to select slices of time, or of frequencies. 
This is a very important service, publicly available, that 
I wish to thank. 
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4. Accretion disks in blazars 


The non-thermal continuum of blazars often dom¬ 
inates in the optical-UV bands, making the accre¬ 
tion disk invisible. And yet, in FSRQs, we do see 
broad emission lines, that should be produced by clouds 
photo-ionized by the disk flux. The old idea of the 
beamed continuum being even stronger in BL Lacs, 
such to hide also the emission lines (besides the disk 
flux) is in general wrong, but it may still be true in a 
few cases. Now we believe that the “genuine” BL Lacs 
intrinsically lack the broad emission lines, because their 
accretion disks are in the low radiative regimes (ion sup¬ 
ported tori 0, ADAF 0, CDAF |70|. ADIOS [11]) 
and therefore have mass accretion rates below a critical 
value in units of Eddington, corresponding to disk lumi¬ 
nosities Ldisk/LEdd$10 -2 [68], [82]. There is then a “di¬ 
vide” in terms of Ldisk/LEdd (or, equivalently, in terms of 
the accretion rate m in Eddington units), distinguishing 
BL Lacs and FSRQ fidll . Since we believe that the cor¬ 
responding parent populations are FRI and FRII radio¬ 
galaxies, they should be characterised by the same di¬ 
vide [36i]. All these issues require the knowledge of 
the mass of the black hole Mbh. There are mainly four 
methods for estimating it. 

re- 


The first is the popular virial method 1 1041] . 
quiring the measurement of the FWHM of the Ha or 
H/j or Mgll or CIV broad lines, and the luminosity of 
the continuum close to the line frequency, if this is not 
contaminated by the beamed radiation from the jet. 

The second methods uses the tight correlation be¬ 
tween Mbh and the velocity di spers ion <x of the galaxy’s 
bulge or spheroid ll28ll . |3JJ. [103]. This is the way in 
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Figure 3. Swift/X RT and NuSTAR, observing simultaneously, cover 
the 0.3-70 keV band. The top panel shows the overall SED of PKS 
2149-306, a blazar at z = 2.345: the optical-UV bump is due to 
the accretion disk, while the extremely strong and hard X-ray emis¬ 
sion is due to the beamed jet emission, and it is interpreted as in¬ 
verse Compton scattering off the photons produced mainly by the 
infrared torus (responsible for the hump at ~ 10 13 Hz). The bot¬ 
tom panel is a zoom on the X-ray portion of the SED, showing 
how the spectrum varied between 2013 December (blue/black points) 
and 2014 April (red/magenta points). Other archival points are also 
shown, to illustrate the variability amplitude of the source. Note that 
in the Fermi/LXY and the source is weak, and the spectrum steep. In 
fact these powerful and high-z blazars emit most of then* luminosity 
around ~1 MeV, and are then better found through hard X-ray sur¬ 
veys rather than >100 MeV surveys (of comparable vF v sensitivities). 
Adapted from f 9^1 . 
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which we can measure Mbh of some relatively nearby 
(z < 0.4) BL Lac objects | 73 |. 

The third method can be applied to relatively nearby 
BL Lacs whose host galaxy light can be distinguished 
and separated from the jet emission [84]. The black hole 
mass can then be found through the relation between the 
bulge luminosity and Mbh [600 ■ 

The fourth method is the oldest and it is based to the 
modelling the thermal continuum with a disk emission 
model. The simplest is a standard, Shakura & Sunyaev 
[85] model, for which the disk emits black-body radia¬ 
tion with increasing temperature for smaller radii B8811 , 
1611. 1621. lf94ll . 11 Kill . The revival of this method T inI 
is due to the realisation that the luminosity of the broad 
emission lines is a good proxy for the disk luminosity, 
that in turn is proportional to the mass accretion rate 
M (through the efficiency 77 defined by Ldisk = pMc 2 ). 
Therefore, in principle, one can find Mbh even when the 
peak of the disc emission is not visible, and also when 
the continuum is partially contaminated by the jet emis¬ 
sion. SwiftlUV OT, with its optical-UV coverage, can 
be used together with WISE or other infrared facilities 
to study the disk emission at its peak (or close to it). 
When the peak is visible the statistical error on the de¬ 
rived black hole mass is small (i.e. less than 50%), in 
comparison to the virial method. 

The top panel of Fig. [4] shows an application of this 
method to the radio-loud source SDSS 0131-0321, at 
z = 5.18 [51|. It shows the IR-UV SED in the rest 
frame and the ocher vertical line indicates the frequency 
of the hydrogen Ly a line. The grey stripe is the estimate 
of the disk luminosity derived by the broad Mgll and 
Ly a lines, taking into account that we observe only the 
non-absorbed part of the latter. The points are archival 
data. Consider that one of the 2MASS data points in¬ 
cludes the Mgll line contribution (and correspondingly 
the point is above the fitting model). 

The solid and dashed lines are three standard disk 
models, with the same Ldisk and slightly different Mbh: 
9, 11 and 14 billions of solar masses. The solid blue line 
shows also the contribution of two two black-bodies at 
different temperatures but similar luminosities, thought 
to be produced by the torus (as found in other, radio¬ 
quiet sources observed by WISE 0 ). 


5. Confirming the “blazarness” of high-z blazar 
candidates 

The simplest way to select sources with jets pointing 
close to our line of sight is to consider sources with a flat 
radio spectrum having high values of radio-loudness 
Rl = E 5 ghz / E 2500 A' where Lsghz and F 2 50oA are the 
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Figure 4. Top panel: the IR-optical data (in the rest frame) of the 
radio-loud quasar SDSS 0131-0321 (z = 5.18), can be well fitted by 
a standard accretion disk model f85ll . plus some contribution from a 
reprocessing infrared torus 0. The grey stripe indicates the range 
of possible vL v peak luminosities of the accretion disk, as derived by 
the broad line luminosities. The black hole mass is (11 ± 2) x 10 9 M Q . 
The small uncertainty is due to the well visible (and constrained by 
the data) peak of the accretion disk component. Bottom panel: entire 
SED of SDSS 0131-0321 with the recent Swift/XKT data. The SED 
can be fitted with a leptonic one-zone model, whose jet is observed 
with a viewing angle between 3 and 5 degrees. Since the viewing 
angle is small, and the jet emission is beamed, there must exist many 
other (hundreds) of similar quasars (with the same black hole mass) 
with jets pointing in other directions From fUli. 
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Figure 5. A sequence of models illustrating the different beaming pat¬ 
tern for the synchrotron and the External Compton (EC) emission that 
occurs in a powerful blazars, following [1251. The very same source 
is observed at different viewing angles 0 V . as labelled. It can be seen 
that the EC luminosity has a stronger dependence on 0 Y than the syn¬ 
chrotron one. The dashed grey lines are the synchrotron component. 
For the chosen T. the proposed divide between blazars and parent pop¬ 
ulation occurs at 0 V = arcsin(l/r) = 4.4°. 

monochromatic fluxes at 5 GHz and at 2500 A, respec¬ 
tively. The rationale of this method is the fact that the 
flat spectrum radio flux is produced by the jet and is 
amplified by beaming, which strongly depends on the 
viewing angle, while the UV flux is not beamed. How¬ 
ever, there are a few caveats: i) it may happen that both 
the radio and the UV are beamed, giving moderate or 
small values of R\\ ii) the UV can be thermal emis¬ 
sion, and therefore quasi-isotropic, but party absorbed 
by dust, giving large /xj also when the viewing angle is 
moderately large. 

Therefore we do not have a clear definition of what 
a blazar is, beyond the generic requirement that the 
jet is observed at “small viewing angles”. But how 
small? We proposed that the natural angular scale divid¬ 
ing blazars from their parent population is the beaming 
angle, namely 1 /T. This definition implies that for each 
observed blazar there are other 2T 2 misaligned sources 
sharing the same intrinsic properties of the blazars, in¬ 
cluding the black hole mass and disk luminosity. This 
also implies that we have to find a good method to ob- 
servationally classify blazar candidates as blazars. Here 
theory comes to help. If most of the high energy emis¬ 
sion is through inverse Compton scattering between rel¬ 
ativistic electrons in the jet and seed photons produced 


outside the jet (External Compton, EC), then the radia¬ 
tion pattern is not isotropic even in the comoving frame. 
This is because in the comoving frame (by aberration) 
most external photons are coming from the forward di¬ 
rection, and head-on scatterings are more energetic than 
tail-on ones. The observer at a small viewing angles 
sees not only the usual beaming pattern, but some ex¬ 
tra emission due to the extra power emitted in the for¬ 
ward direction even in the comoving frame. At large 
viewing angles, instead, the observer will see less radi¬ 
ation. In other words, there is a difference between the 
synchrotron and SSC beamed radiation (usual beaming 
pattern) and the EC emission (beaming+some extra pat¬ 
tern), making the EC flux more sensitive to the viewing 
angle [25]. Fig. 0 shows how the model SED changes 
by changing the viewing angle 6 W . Notice that from 0° to 
6 ° the synchrotron luminosity changes by two orders of 
magnitude, while the EC one changes by three orders. 
Notice also that the X-ray spectrums become steeper 
for larger viewing angles, until the X-ray corona domi¬ 
nates the emission: in the shown example this happens 
for 0 V ^8°. 

This implies that the X-ray flux, relative to the opti¬ 
cal and the radio, becomes a tool to estimate the viewing 
angle of a candidate blazar, if the EC component dom¬ 
inates the observed flux. Fortunately, the EC spectrum 
is usually much harder than the SSC one, and we can 
discriminate. To summarise: if we want to find new 
blazars, we first select flat radio spectra objects with a 
large radio-loudness (i.e. > 100), then we observe it in 
the X-rays. If we see a strong X-ray flux relative to the 
optical and a hard spectrum, then the object is with good 
probability a blazar. The bottom panel of Fig. [4] shows 
a borderline case, because the X-rays, as observed by 
Swift/X RT, are not as strong (relative to the optical) as in 
the PKS 2149-306 (see the top panel of Fig. [3}, but they 
may have a hard spectrum. The source can either be a 
blazar (corresponding to the predicted spectrum shown 
by the green dashed line) or have (9 V slightly larger than 
1 /T (solid blue line). To discriminate we should observe 
the source at higher X-ray energies, i.e. with NuSTAR. 

From what said it should be clear that this method can 
work only with FSRQs, because only in these sources 
we surely have important sources of external photons 
(broad lines and infrared photons for the toms). Fur¬ 
thermore, the emission zone should be located within 
the broad line region, or at a distance from the black 
hole smaller than the toms one, for the EC emission to 
become effective. 
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6. Hard X-ray luminosity function and black hole 

mass function 

Ajello et al. 0] studied all blazars detected by 
Swift/BXT after the first 3 years of operations (see also 
0.0 for more recent BAT catalogs), and constructed 
the corresponding blazar luminosity function. The num¬ 
ber of detected blazars was very limited: 38 in total: 26 
FSRQ and 12 BL Lacs (excluding the Galactic plane 
from the analysis). Of these, 10 blazars (all FSRQs) are 
at z > 2, and 5 at 3 < z < 4. All these 10 blazars 
have a 15-55 keV luminosity L B at > 10 472 0-1 


erg s 


0. We |431 studied these 10 high-z blazars and found 
that the black hole mass of all of them exceeded 10 9 M Q . 
Therefore the number density of blazars at z > 2 with 
Lbat > lO 47 2 erg s -1 can be used to estimate the num¬ 
ber density of heavy black holes at these redshifts. Us¬ 
ing the average value of T we can reconstruct the to¬ 
tal comoving number density of jetted sources having 
Mbh > 10 9 M q . The result shown in Fig. [6]is intriguing. 
The density of active and heavy black holes in jetted 
sources, as derived by using BAT data, peaks at z ~ 4. 
We stress that this is the density of active black holes, 
namely having Ldisk > OTLecM. We can contrast this 
density profile with the one derived by the blazars de¬ 
tected by Fermi/LAT and studied in |5|. shown in Fig. [6] 
as the light blue line. In this case the mass density pro¬ 
file of heavy and active black holes peaks at z ~ 1.5, but 
the peak is at a lower density value. Since we are con¬ 
sidering strongly accreting sources, we expect that most 
of their jet luminosity is emitted in the hard X-rays or 
in the y-ray bands, therefore in the Swift/BAY or in the 
Fermi/LAI bands. We conclude that most heavy black 
holes in jetted sources are indeed formed at z ~ 4. 


By integrating the optical luminosity function of |54] 
above L opt = 10 47 erg s -1 , we can have an estimate of 
the density of heavy and active black holes for radio¬ 
quiet quasars (the luminosity limit corresponds to nearly 
the Eddington luminosity for a 10 9 M G black hole), 
shown by the blue line in Fig. [6] This profile peaks 
at z ~ 2.5 (see Fig. 2 in [[45}] for different L opt thresh¬ 
olds). The conclusion is that very heavy black holes in 
jetted sources form earlier than very heavy black holes 
in radio-quiet sources. Is this because jets are born pref¬ 
erentially in heavy black hole systems, or, on the con¬ 
trary, it is the jet that helps a fast growth of the black 
hole? 


As lolley & Kuncic [57] suggested, the available 
gravitational energy of the infalling matter could be 
used not only to heat the disk, but also to amplify the 
magnetic field necessary to launch the jet. In this case 
the disk is colder, and becomes Eddington limited for 
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Figure 6. Top: Comoving number density of blazars powered by 
“heavy and active” black holes (M > 10 9 M Q , Ldisk/^Edd > 0-1) as 
a function of redshift. The orange stripe is derived by integrating 
the [15-55 keV] luminosity function (LF) (01, as modified in (43h ) 
above Lbat = 2 X 10 47 erg/s, and multiplying the derived density by 
2T 2 = 450 (i.e. T = 15). The light blue stripe is derived by integrating 
the y-ray LF 01 above Ly = 10 48 erg/s. The blue stripe is derived in¬ 
tegrating the LF of radio-quiet quasa rs t54ll above L opt = 10 47 erg/s, 
as labelled (see also f 105h and fl08ll ). The grey stripe, described in 
f43h can be considered as an upper limit. The (yellow) pentagon la¬ 
belled 1023 1146 is the density inferred from the existence of two 
blazars, B2 1023+25 (H, and SDSS J1146+4037 (H] at z > 5 in the 
region of the sky covered by the SDSS+FIRST surveys. Adapted from 
18311 . Bottom: Mbh as a function of time (bottom axis) and redshift 
(top axis). Accretion starts at z = 20 onto a black hole seed of 10 2 M G , 
10 4 M© or 10 6 M o , with different efficiencies, as labeled. The larger 
77d, the smaller the amount of accreted mass needed to produce a given 
luminosity, and the longer the black hole growing time. If part of the 
accretion energy goes into launching a jet, however, < rj and the 
growth time decreases. From l47h . 
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a higher accretion rate. If there is a large reservoir of 
matter that can be accreted, we can have large M and 
relatively small Ljisk even if the total efficiency of the 
accretion is // = 0.3. In practice, the total efficiency 
Tj — //a + //jet is the sum of the radiative efficiency r\& and 
the jet efficiency 7/j et . 

The bottom panel of Fig. [6]shows the black hole mass 
as a function of time (or redshift, upper x-axis) for ac¬ 
cretion limited by the Eddington rate and for 10 2 , 10 4 
and 10 6 M o of the seed black hole, assumed to be at 
z — 20 (the figure is taken from [83]). For the low¬ 
est value of the seed, and for the usual 7/ = 7/d = 0.1 
value, it is not possible to have 10 9 M Q black holes at 
z> 2. Instead we do find larger masses at higher red- 
shifts. Increasing the seed black hole mass does not 
help much: even a 10 6 M o seed reaches 10 9 M Q not ear¬ 
lier than z = 4. We must assume a lower efficiency 7/. 
But this contrasts with the general idea that jets are as¬ 
sociated with black holes spinning rapidly, close to the 
maximum value. In this case the innermost stable orbit 
approaches the gravitational radius R g = GMbh/c 2 , and 
the efficiency is close to the maximum value (which is 
77 = 0.42 in principle, but only 7/ = 0.3 in reality, see 
1102]). The problem can be solved if we assume that 
even in 7/ = 0.3 systems, most of the gravitational en¬ 
ergy is used to amplify the magnetic field necessary to 
tap the rotational energy of the hole. The “disk effi¬ 
ciency” 7/d can then be smaller than 0.1: to produce a 
given disk luminosity, the mass accretion rate is cor¬ 
respondingly larger, and the black hole can then grow 
faster (blue dashed lines in Fig. [6|. 

But why the number density of heavy black holes 
derived from the y-ray luminosity function peaks at a 
much smaller redshift, similar to the peak of the radio¬ 
quiet ones? Are these black holes in their growth age, 
or are they resurrected black holes, e.g. by a merging 
event? We do not know yet. 

In any case, there seems to be two epochs of forma¬ 
tion of heavy black holes. One, at z ~ 4, in which the 
majority of Mbh > 10 y M o with jets is born and grows, 
the other, at z ~ 2 where the majority of heavy black 
holes in radio-quiet quasars are born. Since we are con¬ 
cerned only with heavy black holes, that are a minority 
of the total, this might not affect the radio-loud frac¬ 
tion of the entire population of quasars, that is seen to 
decrease with redshift (but not with luminosity) [53]. 

These results are a direct consequence of the 
Swift/ BAT survey. The fact that it includes blazars at 
higher redshifts than FermifLXT is the consequence of 
two effects: i) increasing the redshift, the Compton peak 
of the SED moves closer to the observed high energy X- 
ray band, and ii) more powerful sources have the Comp- 



Figure 7. Cartoon of the possible location of the emitting region. 
The most economic location is within the BLR (1): there the inverse 
Compton process can benefit of the photons produced by the BLR. If 
the location of the emitting region is beyond the BLR (2), but within 
^torus> there still is the benefit of using the infrared photons produced 
by the torus. Beyond R t oi- U s (3) the energy density of the external ra¬ 
diation drops rapidly, and the main source of photons for scattering 
is the internally produced synchrotron radiation. The emitting region, 
in this case, must be a small part of the jet, otherwise the variability 
timescales become too long. 


ton peak located (in the rest frame) at somewhat smaller 
energies (i.e. close to 1 MeV), in general agreement 
with the blazar sequence (see Fig. 11 of 14311). 


7. Location of the emitting region 

There is an ongoing important debate about the loca¬ 
tion of the emitting region, and if the one-zone model 
is an over-simplification or instead if it is a good repre¬ 
sentation of reality. 

First, we have to be aware that the jet has surely 
several emitting regions: we do see some of them in 
the VLBI observations, and it is very likely that this 
“knotty” structure exists also at the smallest scale. On 
the other hand, if the jet is accelerating, the beaming an¬ 
gle of the most inner parts is very large, implying that 
the corresponding luminosity is diluted in a large solid 
angle, and therefore not_visible in blazars (but it could 
be in radio-galaxies, |4^]). However, there could be 
several emitting zones also when F has reached its final 
value. After all, this is what predicted by the “internal 
shock model” for blazars [791], 19ill . And yet we often 
see coordinated variability at different frequencies, im¬ 
plying that there is one region that dominates the bolo- 
metric luminosity. Swift gave a fundamental contribu¬ 
tion to this issue. 
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Within the BLR — The variability timescales f var 
are often short, of the order of few hours [98], and 
imply rather small dimensions of the emitting region, 
R<ct val 6/(l +z)~4x 10 15 (f var /3h)(<5/10)/(l + z) cm 
(5 is the relativistic Doppler factor). In turn, for conical 
jets of aperture angle if/ ~ 0.1 rad, we derive that the 
emitting region is at /?diss ~ 4 x 16 cm from the black 
hole. Since the broad line region is supposed to be at 
a distance Rblr ~ 10 17 (Ldisk/10 45 erg/s) 1 ^ 2 cm, most 
FSRQs should have their emitting region preferentially 
within the BLR (zone 1 in Fig. |7J. In this case there 
should be a well defined signature in the y —ray spec¬ 
trum of blazars, since the photons of the broad lines at 
high frequency (e.g. Hell) could be the targets for the 
y—y —» e 1 process. This should produce an absorption 


feature in the high energy (GeV) spectrum I174ll . 119311. 

Beyond the BLR but within the torus — Another 
possible location for the most active part of the jet of 
FSRQs is outside the BLR, but within R toras [12]. By 
“torus” we mean a reprocessing region similar to the one 
in radio-quiet quasars, intercepting 10-40% of Ldisk and 
re-emitting it in the IR, probably with a multi tempera¬ 
ture blackbody spectrum, the hottest being T ~ 2000 K, 
the sublimation temperature of the dust. This is the “2” 
zone in Fig. [7] Both Rblr and R t orus should scale as 
^kiisk ■ This i m pli es that within Rblr the radiation en¬ 
ergy density Ublr is constant. Also between Rblr and 
R torus the radiation energy density f/ tor us is constant but 
with a reduced value: 6 'bi.r > U toms . Since it is likely 
that the magnetic field B scales as R 1 , corresponding to 
a constant Poynting flux (oc R 1 \ ' 1 B 1 ), an emitting zone 
further out in the jet implies a less magnetised zone. A 
given Compton to synchrotron luminosity ratio can usu¬ 
ally be produced in two well defined zones, one within 
Rblr and one between Rhi.r and R 1(irus [Bj/]. To choose, 
apart from the expected variability timescale, we can 
use the peak frequency of the Compton hump: if smaller 
than a few MeV it may indicate the need for IR seed 
photons. 

Beyond the BLR and the torus — In some blazars 
there seems to be a correlation between y-ray flares, ra¬ 
dio flares, and a switch in the polarization angle [58], 
[64]. This suggests that the emitting region is located 
parsecs away from the black hole, in the VLBI re¬ 
gion. This contrasts with the short f var , and also with 
the cooling timescale, that should be long at that dis¬ 
tances (small magnetic field, no important sources of 
external photons). And yet, the ultrafast variability seen 
in the TeV band is a puzzle. Two TeV blazars, Mkn 
501 and PKS 2155-304 showed significant (factor 2) 
variations of their TeV flux in f var =3-5 minutes [1], 
llBll . This is already difficult to explain 0], fBBIl . IMF ], 


[ 65], but the real puzzle came when PKS 1222+216, at 
z = 0.431, was also observed to vary, at a few hundreds 
GeV, in f var ~ 10 minutes (2l, corresponding to a size 
R < c/ val <5/(1 +z) - 5 x 10 14 (5/20) cm. PKS 1222+216 
is a FSRQ with broad emission lines, and this extremely 
small TeV emitting region cannot be located within the 
BLR, whose photons would absorb the emission above 
~20 GeV. It must be located outside. This remains true 
even if the BLR has a flattened geometry (as suggested 
by e.g. & IBbll and ll23ll ). The TeV emitting region of 
PKS 1222+216 can be located between Rblr and R torus , 
but must be much smaller than the cross sectional ra¬ 
dius of the jet at these distances, or the jet itself must 
shrink, due to strong recollimation and focusing of the 
flow (e.g. [92]; [16], [67]). Alternatively we may have 
complex reconnection events [[53]; or even a photon to 
axion transition, to survive the y—y process [99]. 


8. Jets and accretion 


In their pioneering work [78] Rawlings & Saunders 
found that the radio-lobes of AGNs, to exist, require 
an average power that is of the same order of the ac¬ 
cretion disk luminosity. Uncertainties were large, since 
there was no idea of the contribution of the protons to 
the total lobe energetics. Attempts to measure the jet 
power continued, modelling the SED using the VLBI 
radio size and a limit on the 5-factor from the require¬ 
ment not to exceed, by the SSC process, the observed 
X-ray flux lUl. [201. 

In a similar way, modelling the radio-optical and X- 
ray emission at relatively large jet scale, as resolved 
by Chandra, allowed to estimate the jet power at these 
scales [37], and to compare it with the values found for 
the compact jet [97]. 

Another more recent way to indirectly measure the 
jet power is through X-ray cavities seen in (relatively 
nearby) radio-galaxies, coinciding with their radio¬ 
lobes: The mechanical power P cav of the cavity can be 
derived by setting P cav = PV/t age , where f age is the rise¬ 
time of the cavity, P is the pressure and V is the volume 
of the cavity [27]. 

With the advent of CGRO and now Fermi we at last 
knew that most of the power of blazars was emitted at 
high energies. In powerful FSRQs, the inverse Comp¬ 
ton luminosity can be larger than the synchrotron one 
by up to 3 orders of magnitude. We can then derive 
the power spent by the jet to produce the observed jet 
bolometric luminosity Lj et (corresponding to L' et in the 
comoving frame). If the processes are synchrotron and 
SSC (isotropic in the comoving frame) we have, for one 
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jet: 


where the last equality is for viewing angles 0 V ~ 1/r, 
for which d = F. This is the entire power carried by 
the produced radiation, at all angles. If most of the lu¬ 
minosity is produced by the external Compton process, 
the beaming pattern is not oc 5 4 , but oc d 4 (d/T) 2 [125], 
lf32h . Setting (L' et ) as the angle-averaged luminosity in 
the comoving frame, we have/or one jet |44|: 


Pi ad 


<^et> 

4n 


f 

J An 


5 6 {0) 


da 




jet/ 


1 6 T 4 Lje, 1 6 Ljel 

5 6 6 {0m) 5T 2 


( 2 ) 


Again, the last equality is valid for 6 V ~ 1/T. We see 
that the difference between the two cases, for blazars 
(i.e. when 6 ~ T), is only in the numerical coefficient. 

To estimate the power P ra( j we need to know the bolo- 
metric jet luminosity and the value of T. The latter can 
be derived by modelling, or by assuming that the values 
derived from the apparent superluminal velocity occur¬ 
ring at the VLBI scale (i.e. 10 pc) are the same of the 
emitting region. It is therefore a very robust quantity, 
almost model-independent. On the other hand, P rM j is 
only a lower limit to the real jet power: if all of the jet 
power is spent to produce the radiation we see, then the 
jet should stop, and could not produce the superluminal 
blobs or energise the extended radio structure. In low 
powerful, TeV BL Lacs, something of this kind may in¬ 
deed happen [38]: these are sources requiring the high¬ 
est T in the emitting region, and yet they do not show 
superluminal motion, nor strong extended structure. 

We want to compare P m d with the disk luminosity 
Tdisk- We then need a sample of blazars detected in the 
y-ray band ( where most of Lj et is emitted) and for which 
we can reliably estimate L ( |,, k . This sample has been as¬ 
sembled rather recently by Shaw et al. [86;], [87], that 
observed spectroscopically hundreds of blazars (both 
FSRQs and BL Lacs) detected by FermifLACV. 

We H, @ have selected all the 217 objects show¬ 
ing at least one broad emission line in their spectrum. 
Among those we have 26 blazars classified as BL Lacs 
according to the classical definition (i.e. equivalent 
width smaller than 5 A) but that indeed have a broad 
emission line in their spectrum. Therefore they belong 
to the low line luminosity tail of FSRQs rather than be¬ 
ing classical BL Lacs. We constructed the overall SEDs 
of all objects, calculated the bolometric jet luminosity 
Lj et , and applied a one-zone model to fit the data from 



L d isk [ er g s_1 ] 

Figure 8. The power spent by the blazar jet to produce the radi¬ 
ation we see, P ra d is comparable with the accretion disk luminos¬ 
ity Tdisk- All shown blazars have broad lines, used as a proxy for 
the disk luminosity, and are detected by Fermi/LXT, and thus have a 
well determined bolometric non-thermal luminosity L , from which 
P r ad ~ Ljet/r 2 is derived. Therefore the result P ra d ~ Ldisk is quasi 
model-independent, since T is rather well determined. The plotted 
P ra d refer to two jets. Adapted from |f48l] , with the addition of all 
known z > 4 blazars. 
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the mm to the high energy y-rays. We could then de¬ 
rive the bulk Lorentz factor F for each source, finding a 
narrow distribution (10 < T < 17). This agrees with the 
values found through superluminal motion. Therefore 
we could reliably estimate P t a< j for all sources: at this 
stage, the applied model is used only for estimating T. 

The disk luminosity could be estimated through the 
observed emission lines, using the average templates of 
0 °r Goi and then multiplying by the average ratio 
Ldisk/F blr ~ 10, or using the value directly listed in 
I8j|], or using the value found fitting the IR-UV SED 
with a Shakura & Sunyaev disk spectrum. We checked 
that the three methods agreed. In the few cases for 
which they did not, we choose the disk-fitting value. 
Fig. [8] shows the result: P a d ~ Ldisk with a probability 

P < 10. 8 to be random, even taking into account the 

common redshift dependence. The coloured stripes in¬ 
dicate 1, 2, and 3cr (vertical) dispersion (cr =0.5 dex). 
Just for curiosity, we plot also the values for a sample 
of z > 4 FSRQs having good optical-UV spectral cov¬ 
erage (hence a reliable Ldisk) and X-ray data, used to 
estimate P m <±. Since they have not (yet) been detected 
by Fermi/LXT, their inverse Compton hump luminosity 
is uncertain, and their real Lj et could be larger. From the 
found correlation we can conclude that the lower limit 
on the jet power is of the same order of the disk lumi¬ 
nosity, and therefore that the real jet power should be 
larger than Ldisk- 

The problem is to find the real Pj et . In the past, there 
have been four main uncertainties concerning the esti¬ 
mate of P j et by using emission models to fit the spectra: 


1. The emission model could be hadronic, rather that 
leptonic. 

In this case one needs very energetic protons, that 
pair-produce in photo-meson interactions. The jet 
power, as calculated in 115] is greater than in lep¬ 
tonic models. 


2. The emitting region could be far away, with SSC as 
the most important emitting process. 

Usually, a simple SSC does not fit well the SED 
of FSRQs. Besides, if the emitting region is be¬ 
yond Rtorusi there are no important sources of ex¬ 
ternal photons, and presumably the magnetic field 
is smaller. This means that the synchrotron and in¬ 
verse Compton scattering process become less ef¬ 
ficient radiators. To produce the radiation we see, 
we need more particles, and the jet kinetic power 
is larger. 

3. The assumption of one proton per electron makes 
the proton kinetic energy dominating Pj et . But most 
of the electrons, for energy distribution N(y) oc y~P 


with p > 1, lies at the low energy end. If there is 
low energy cut-off, we might not notice it in the 
data. 

In FSRQs, the EC process dominates, making the 
soft X-ray spectrum very hard. This part of the 
SED is produced by low energy electrons scatter¬ 
ing either broad line photons, mainly of the hydro¬ 
gen Lya line, or IR photons from the torus. We do 
have control of the low energy tail of the particle 
distribution. 

4. The assumption of one proton per electron could 
be wrong because of electron-positron pairs. 

Pairs cannot be produced (in appreciable number) 
in the emitting region, otherwise they would repro¬ 
cess the spectrum, especially at the low X-ray fre¬ 
quencies, making the X-ray spectrum softer than 
observed. There is the possibility to produce pairs 
at the base of the jet, in a region with small T, 
whose radiation is overwhelmed by the main emit¬ 
ting region in blazars, but could be visible in ra¬ 
dio galaxies. However, the SED is required to be 
finely tuned, to produced a sufficient number of 
pairs 0, and this is unlikely. See also 0, JH, 
|44] for additional arguments concerning the pres¬ 
ence of pairs. 

In addition to these arguments, and in completely in¬ 
dependently way we also have the result of Nemmen 
et al. ifnll . that estimated P ra d and Pje\ for blazars and 
Gamma Ray Bursts, finding that they lie on the same 
correlation, indicating Pj et ~ 1 0P rdi \. We recover their 
result if we assume one proton per electron. 

Fig. E shows how the different forms of powers are 
distributed: If, is the Poynting flux. If is the kinetic 
power of the emitting electrons (thus including their rel¬ 
ativistic random energy), and P p is the kinetic power of 
protons, assumed to be cold in the comoving frame. One 
can see that to produce P a &, the Poynting flux and the 
electron kinetic power are not sufficient. One needs an¬ 
other reservoir of power, and the simple assumption is 
that this is provided by protons. 

Finally, Fig. ITOlshows the total jet power as a function 
of M, The latter is found assuming that Ljj sk = r/Mc 1 2 3 , 
with p = 0.3, i.e. for a maximally efficient accretion and 
maximally rotating black hole. The yellow stripe indi¬ 
cates equality, while the black line is the best fit of the 
correlation. We find that Pj e t>Mc 2 , and yet it is corre¬ 
lated with it. This is an apparent paradox: the fact that 
it correlates let us think that the jet is powered by ac¬ 
cretion, while the fact that its power is greater than Me 2 
implies that this is not possible. 

The solution of this paradox is the following: part of 
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42 44 46 48 


Log Power [erg s -1 ] 

Figure 9. Distribution of the different forms of jet powers compared 
to the disk luminosity. The jet powers assumed two jets. The darker 
(blues) hatched regions correspond to the “BL Lacs” of the sample. 
P p and P e are the kinetic power of emitting electrons and of the (cold) 
prot ons, Pq is the Poyinting flux, Ldisk is the disk luminosity. From 



Me 2 [77 = 0.3, erg s _1 ] 

Figure 10. The total jet power Pj et as a function of the mass 
accretion rate, assuming an accretion efficiency 77 (defined through 
Ldisk = pMc 2 ) e qual to 0.3, a value appropriate for maximally spin¬ 
ning black holes Il02ll . The yellow line is the equality line, while the 
black line is the best fit. The jet power is of the same order of and 
possibly larger than, Me 2 . Therefore, despite the strong correlation 
between the jet power and the accretion rate, the latter is not enough 
to power the jet. Another source of power is needed, such as the ex¬ 
traction of the rotational energy of the spinning black hole. Adapted 
from jf48h . with the addition of all known z > 4 blazars. 
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the gravitational energy of the infalling matter is used 
to amplify seed magnetic fields up to equipartition with 
the mass energy density ~ pc 2 of the matter accreting 
at the rate M. We then have B 2 oc p oc M. According 
to the Blandford & Znajek process jsj], the jet power 
depends on ( aMB ) 2 , the square of the product of the 
hole spin, mass and the magnetic held at the horizon. 
Therefore Pj et oc a 1 B 1 oc crp oc M. This explains why 
P j e t correlates with the accretion power. 

To explain why Pj et >Mc 2 we are forced to assume 
that the jet production process extracts energy not only 
from accretion, but mostly from the rotational energy 
of the black hole. This is done by the magnetic held 
that must be thought as a catalyst for the process. Since 
Pj e t>Mc 2 > L ( ii. s k, the extraction of the rotational en¬ 
ergy must be extremely efficient. This is fully consistent 
with general relativistic magnetohydrodynamic numer¬ 
ical simulations EH in which the average outflow¬ 
ing power in jets and winds reaches 140% of Me 2 for 
dimensionless spin values a =0.99. Occasionally, the 
magnetic energy density can exceed the energy density 
~ pc 2 of in the vicinity of the l ast s table orbit, and the 
accretion is temporarily halted EH EH EH pro¬ 
viding a way to explain the observed variability. 
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